Phase transition in exotic nuclei along the N = Z line 
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The abrupt structure change from the nuclei of A'^ = Z < 35 to those of A'^ = Z > 36 is investigated 
by means of shell model calculations. The basic features of the even-even and odd-odd nuclei under 
consideration are nicely reproduced. A sudden jump of nucleons into the upper 39/2^5/2 shell at 
A'^ = Z = 36 is found to be the main reason that causes the qualitative structure difference. It is 
argued that the structure change can be viewed as a decisive change of the mean field, or a phase 
transition, along the N = Z line. 

PACS numbers: 21.10.-k, 21.10.Re, 21.60.Cs 



Phase transition in a many-body system refers to an 
abrupt, qualitative change in wave function or mean 
field. This subject is of common interest for many sub- 
fields. Nuclei, being finite quantum systems composed of 
strongly correlated protons and neutrons, uniquely dis- 
play transitional features. Often discussed in nuclei are 
two types of phase transition: the phase transition to 
superfluidity and to deformation. The first one is easy 
to trace with the standard BCS theory which gives the 
critical condition The second one, the transition to 
deformation, takes place as some control parameters vary 
along an isotopic (isotonic) chain, which, for example, 
brings the system from a spherical to a deformed region. 

The second type of phase transition, also called shape 
phase transition has been theoretically studied but 
mainly by means of algebraic models @. Well-suited 
cases are those described by an algebraic Hamiltonian 
with a dynamical group, where a transition-driving con- 
trol parameter appears explicitly in the Hamiltonian. 
Large-scale shell model descriptions 4, 5] for deformation 
in medium-mass nuclei have become possible in recent 
years. The shell model can provide a more fundamental 
basis to the study of shape phase transition. In fact, the 
Monte Carlo Shell Model calculation Q has shown the 
first example. The advantage of a shell model study is 
that one may see microscopically the origin of a phase 
transition by analyzing the wave functions if transition- 
driving mechanism is contained in the effective interac- 
tion. The aim of the present Letter is to carry out such 
a study using the spherical shell model, to show a shape 
phase transition in A ~ 70 nuclei along the N = Z line, 
and to understand the cause of the transition by studying 
the occupations in the nucleon orbits. 

The N = Z nuclei around ^ 70 exhibit several 
unique phenomena 0, Q, and therefore have attracted 
many theoretical and experimental studies. This is pos- 
sibly due to the rather strong proton-neutron correla- 
tions in these nuclei since protons and neutrons occupy 
the same orbits. In particular, there is apparently an 
abrupt change in structure when the proton and neu- 



tron numbers cross N — Z — 35, which can be clearly 
seen from the graph of spin J versus angular frequency 
uj = {E{J)-E{J-2))/2 (the so-called J-uj graph). Fig. 
[1] shows a separation of two groups for the N = Z nuclei 
from ^^Ge to ^^Sr, one with larger w's corresponding to 
smaller moments of inertia and another with smaller w's 
corresponding to larger moments of inertia. In the left 
group, ^^Sr takes a straight line whose extension inter- 
sects the origin of the J — lo graph. Thus ^^Sr behaves 
like a rotor with a large and approximately constant mo- 
ment of inertia J juj. The J — uj curves for ^^Kr and ''^Rb 
show that these nuclei have a structure rather similar to 
that of ^^Sr although the lowest state of ^^Kr appears to 
be peculiar (see discussions below). In contrast, the nu- 
clei in the right group with N ~ Z < 35, which resemble 
each other, have a moment of inertia only approximately 
half of the values of the left group. This sharp differ- 
ence suggests a qualitative structure change between the 
nuclei with iV = Z < 35 and iV = Z > 36. No nuclei 
sit in between, implying that the structure change with 
nucleon number is sudden. 
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FIG. 1: (Color online) The known experimental data shown 
in the graph of spin J versus angular frequency w. 
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The above structure change along the N = Z Hne rnay 
be comparable to that near = 40 in Ge isotopes [1], 
which has been observed in the (p, t) cross sections and 
B{E2]Q'l 2f) values. For example, ^^Ge has an un- 
usually low second 0+ state below the 2^ state, and so 
does ^^Kr [lo|. In a separate study for Ge isotopes [ll|, 
we have discovered the possible sources for the structure 
change near = 40. The study suggests that the change 
is caused mainly by a sudden jump of nucleons into the 
(7g/2 orbit. We may expect that the gg/2 occupation is 
also the leading source for the structure change along 
the iV = Z line. 

In the early study [ll|, we found it difficult to get 
a sufficient (79/2 contribution within the p/5/239/2 shell 
model, in spite of the fact that this model s pac e is capa- 
ble of describing the nuclei ^''Ge and ^*Se ^M,- On the 
other hand, the Shell Model Monte Carlo calculations [l3| 
suggest that the orbit has a cooperative effect that 
enhances the 179/2 contribution. Thus the pfz/29'j/2dh/2 
shell model can be a hopeful model to investigate the 
structure change shown in Fig. [TJ This model space is, 
unfortunately, too large to perform a shell model calcu- 
lation. The model space that we can presently handle 
is a truncated one /5/2?'i/259/2C^5/2- (Note that calcula- 
tion for ^^Sr can be performed by using the extrapolation 
method [l3|-) Setting nucleons in the P3/2 orbit inactive 
is a severe restriction since correlations in the fp shell 
are strong. Nevertheless, this can be compromised by 
adjusting effective interactions. The choice of the present 
model space seems to have grasped the basic physics, as 
we discuss below. 

With the extended P + QQ Hamiltonian [ll|, , we 
have searched for suitable parameters that can reproduce 
the experimental data in this mass region. The level 
schemes shown in Fig. [5] are obtained by using the fol- 
lowing set of parameters. Single-particle energies for the 
/5/2J Pi/2j 59/2i and orbits are 0.0, 0.3, 1.5 and 2.0 
MeV, respectively. For interaction strengths, we take 

go = 0.25(68M), 32 = 0.16(68/^)^/3^ 

X2 = 0.14(68/^)^/3, X3 = 0.04(68/^)2 in MeV. (1) 

The average T — monopole field and T — 1 monopole 
corrections (in MeV) are fixed as follows: 

k° = Hl=^{a,h) = -0.63(68/A) for arbitrary (a, 5), 



(/5/2,Pl/2) = -0.6, 



T=li 



T=l, 



(/5/2,.99/2) - ^^mc (/5/2, 6^5/2) = "0.5, 

(59/2,59/2) = 0.2,i/l=i(d5/2,d5/2) = -0.18.(2) 



It is found that to excite nucleons more efficiently 
to the upper orbits (59/2, '^5/2)7 we should take a 
small ^5/2 single-particle energy close to the gg/2 one 
and large monopole corrections H^'^^{f^f2, 99/2) and 
^^mcr^(/5/2j c?5/2). The above parameter set may be 
named "the effective interaction of the small model" . 

Calculated level schemes for the N = Z nuclei ^*Se, 
^^Kr, ""^Sr, and the N = Z + 2 nucleus ''"'Kr are compared 
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FIG. 2: (Color online) Comparison of experimental and cal- 
culated level schemes for ®*Se, ^^Kr, ^'*Kr, and ^®Sr. 



with experimental data in Fig. [21 Both the ground band 
(in black) and the first excited 0+ band (in blue) are 
shown for each nucleus. The truncated space without the 
P3/2 orbit cannot perfectly reproduce the data. However, 
the qualitative feature for each nucleus is well described. 
The ^^Se nucleus is seen to have a different structure with 
much larger energy intervals than the other three nuclei. 
What we find remarkable is that the calculation correctly 
gives the first excited Oj state in ""^Kr below the 2f one, 
in agreement with the experimental finding [iq| . 

To trace the microscopic origin of the structure 
changes, we calculate nucleon occupation numbers (ua) 
in respective orbits. As we employ the isospin invariant 
Hamiltonian, the proton and neutron occupation num- 
bers in iV = Z nuclei are equal to each other, i.e., 
("a) — ("a)- To simplify the notation, we abbreviate the 
lower orbits (/5/2,Pi/2) and the upper orbits (59/2,^^5/2) 
to "/p" shell and "gd" shell, respectively. We suppose 
that the former represents the fp shell and the latter the 
gd shell. 

The calculated occupation numbers are shown in Fig. 
3. As seen in Fig. 3a, nucleons in ^*Se have a negligible 
occupation in the "gd" shell. Although the Oj state has 
a different configuration from the states of the ground 
band, all of the states are basically constructed in the 
"/p" shell. We should mention that the description for 
^*Se with the present model space may be too simpli- 
fied. However, a more realistic calculation for ^^Se in 
the pf 5/299/2 shell (including the P3/2 orbit) [l^ does 
not give a structure much different from the present one. 
We have confirmed this by improving the p/5/259/2 shell 
model so as to reproduce well the observed coexisting 
oblate and prolate bands [lH|- We thus conclude with 
confidence that the ^^Se nucleus does not essentially oc- 
cupy the "yd" shell. 

In contrast to ®*Se, approximately two protons and 
two neutrons jump into the "gd" shell in the ground state 
Oi and the first excited state 2^ in ^^Sr (see Fig. 3c). 



3 




FIG. 3: (Color online) Proton and neutron occupation numbers in ''^Se, ^^Kr, and '^^Sr. 



Similarly, nucleons in Kr also start to occupy the "gd" 
shell (see Fig. 3b). These results indicate that for ^^Sr 
and ^^Kr, the 0^ state instead of the ground state has 
the ordinary type of configuration, i.e., the occupation 
number deceases gradually as the single-particle energy 
increases. The ground state has the dominant configura- 
tion ("/p")"^''("(;(i")"', the probability of which is about 
0.45 in ^^Kr and is more than 0.7 in ^^Sr. 

The wave functions of the 0^ and Oj states in ^^Kr 
have the following weights for the leading configurations: 

|0+) : 0.22("/p")' + 0.45("/p")'('V")^ + ---, 
|0+) : 0.47("/p")* + 0.40('7p")'(V")' + ---- (3) 

These numbers show that the ground state is constructed 
starting from the excited configuration ("/p")''("5d")'* 
and the second 0+ state is constructed starting from 
("/p")*- The two 0+ states must reverse their order in 
energy if interactions were increased gradually in a treat- 
ment of the perturbation theory. This reversal manifests 
itself in the occupation numbers in Fig. 3b; the 0^ state 
has less nucleons in the "fp" shell and more nucleons in 
the "f/d" shell as compared to the 0^ state. As shown 
in Fig. 3c, ^^Sr has the same characteristic occupation 
numbers obtained with the extrapolation method. 

The calculation thus reveals a large difference in occu- 
pation numbers between ^*Se and ^^Kr, which correlates 
with their qualitatively different picture in the J — w 
graph in Fig. [1] It has been known that states corre- 
sponding to both oblate and prolate shapes coexist in 
68Se UJ. Now our result that the ground state is con- 
structed mainly in the fp shell suggests that the mean 
field of ^*Se is deformed but does not lie very far from 
the spherical one of shell model. In contrast, the ^^Kr 
ground state cannot be described in the perturbation the- 
ory based on the "/p" configuration. This means that a 
decisive, sudden breaking of the spherical mean field of 
shell model occurs when going from ^^Se to ^^Kr. We 
may call the abrupt change in mean field "phase transi- 
tion" . Then what mean field is formed after the breaking 
of the spherical shell model one? The occupation of the 
"gd" shell in ^^Kr and ^^Sr implies a definite formation 



of the Nilsson orbits in the deformed mean field. The 
transition to ''^Kr must be a shape phase transition. In 
fact, the level scheme of ^®Sr shows a rigid rotor behavior 
as seen in Fig. [1] which is qualitatively different from the 
nuclei with N — Z <i^. It has been discussed with the 
deformed mean field language 0, [3 that there is a rapid 
change in deformation with the deformation parameter /3 
from ~ 0.25 in ^SSe, ~ 0.4 in ^^Kr, and ~ 0.5 in ^^Sr. 
The gg/2 orbit has been considered to drive deformation 
in this mass region 

It has been found that in ^^Kr, the nearby 0^ and Oj 
states coexist with different shapes: the excited one is 
considered as a shape isomer [l^, [iMl ■ In our wave func- 
tions ([2]), the 0^ state has a very large component of the 
{"fp")^{"gd")'^ configuration which is the main compo- 
nent of the 0^ state, while the 0^ state has a large com- 
ponent of the ("/p")* configuration which is the main 
component of the 0^ state. The situation that these two 
0"*^ states mutually have both configurations with con- 
siderable amount in the wave functions suggests a rather 
strong coupling between the two states. Such a strong 
coupling can explain why the moment of inertia of the 
2^ state in ^^Kr deviates from the regular rotor behavior 
(see Fig. I); it is simply because the 0^ energy is pushed 
down by the coupling. 

To test the model's reliability, we should also calculate 
the odd-odd nuclei ''"Br and ''^Rb. Calculated energy 
levels are compared with experimental data [l9l . [20| in 
Fig. m The agreement between calculation and experi- 
ment is satisfactory for the restricted model. It is found 
that the configuration ("/p")® is dominant in the ground 
state of ™Br. Contrary to this, for ^''Rb the configu- 
ration ("/p")^("g(i")'* is leading in the low-lying collec- 
tive states in the T = I band (0+,2+,4+) and T = 
band (3+,5+,7+). These results are consistent with the 
clear difference in their moments of inertia between ^"Br 
and '^'^Rb (Fig. [T]). We thus reach the same conclusion 
that the structure of ^''Rb is qualitatively different from 
that of ^°Br and a shape phase transition takes place at 
Z = N = 35. 

The structure change along the Ge isotopic chain man- 
ifests itself also in a rapid increase of B{E2; 0^ — > 2^) 
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FIG. 4: (Color online) Experimental and calculated level 
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FIG. 5: Calculated B(S2)T values for A = 68 - 76 nuclei. 

[Tll |. Similar increase of i3(£'2)t is expected when iV ex- 
ceeds 35 along the N = Z line. In Fig. [5] the calculated 
B{E2)'\ values from ^^Se to ^^Sr clearly show such a fea- 
ture. Here, we used the effective charges Ctt = 1.5e and 
= 0.5e. From the above analysis, the rapid increase 
in B{E2)'\ is attributed to the fact that more degrees of 
freedom are opened up for the quadrupole correlation as 



nucleons occupy the "gd" shell. Thus the present results 
support our previous prediction that the notable in- 
crease of B{E2)^ at iV = 40 in Ge isotopes is caused by 
a jump of nucleons into the "gd" shell. 

We note that our truncated model space exposes its 
weakness in reproduction of absolute value of B{E2)'\. 
For '^^Kr, the calculated value 0.17 e^b^ is only about 
one third of the observed one ^ 0.5 e^6^ [l^. However, 
we know the possible reason for the discrepancy. The 
shell model for ^*Se with the pf 5/299/2 space [1^ gave 
the B{E2)1 value 0.16 e^b^, which is about three times as 
large as the present result 0.054 e^b^. Including the P3/2 
orbit, which contributes to correlations in the fp shell, 
will enhance the B{E2)1 values considerably. For ^^Se, 
while the pf 5/299/2 shell model [ill reproduced a positive 
spectroscopic quadrupole moment {Qsp) corresponding 
to an oblate shape for the 2f state, the present model 
fails. This suggests that the lack of the orbit is 
a serious problem for those quantities sensitive to the 
details of wave functions. Obviously it requires a larger 
model space to describe these quantities. 

In conclusion, the experimental level schemes of the 
A ^ 7Q N ^ Z nuclei show clearly an abrupt struc- 
ture change at TV = Z = 36. We have investigated this 
problem by means of large-scale shell model calculations. 
In spite of the limit of the practicable model space, we 
have been able to reproduce the basic features of the 
level schemes. Through analyzing the wave functions, 
we have concluded that the qualitative structure differ- 
ence between the nuclei with = Z < 35 and those with 
= Z > 36 has the origin of the upper "gc?" shell occu- 
pation. It has been discussed that the structure change is 
caused by a decisive breaking of the spherical shell model 
mean field and a formation of deformed mean field. In 
this sense, we have witnessed a basic mechanism of the 
phase transition to deformation in nuclei. We note that 
calculations in a space larger than pf 5/299/2d5/2 must be 
carried out to describe deformation properly. 
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